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Abstract 

The purpose of this paper is to enhance a correspondence between the dynamics 
of the differential equations y{t) = g{y{t)) on R and those of the parabohc equations 
ii = Aii + /(x, n, Vn) on a bounded domain 0. We give details on the similarities of 
these dynamics in the cases d = 1, d = 2 and d > 2> and in the corresponding cases 
f2 = (0, 1), = T^ and dim(r2)> 2 respectively. In addition to the beauty of such a 
correspondence, this could serve as a guideline for future research on the dynamics of 
parabolic equations. 
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1 Introduction 

In this paper, we want to point out the similarities between the dynamics of vector fields 
in W^ and those of reaction-diffusion equations on bounded domains. More precisely, we 
consider the following classes of equations. 

Class of vector fields 

Let d > 1 and r > 1 and let g G C"(M°',M°') be a given vector field. We consider the 
ordinary differential equation 

y{t) = g{y{t)) t>0 

y(0) = 1/0 G M'^ ^ ■ ^ 

where y{t) denotes the time-derivative of y{t). 

The equation (11.11) defines a local dynamical system Tg{t) on M.'^ by setting Tg{t)yQ = y{t). 

We assume that there exists M > large enough such that 

Vy G M^ \\y\\ >M ^ {y\giy)) < . 

This condition ensures that Tg{t) is a global dynamical system. Moreover, the ball -8(0, M) 
attracts the bounded sets of M*^. Therefore, Tg(t) admits a compact global attractor^ Ag. 
The attractor Ag contains the most interesting trajectories such as periodic, homoclinic 



^To make the reading of this article easier for the reader, who is not familiar with dynamical systems 
theory or with the study of PDEs, we add a short glossary at the end of the paper. 



and heteroclinic orbits^ and any a— or w— limit set^. Therefore, if one neglects the tran- 
sient dynamics, the dynamics on Ag is a good representation of the whole djTiamics oiTg[t). 

Class of scalar parabolic equations 

Let d' > 1 and let Q be either a regular bounded domain of M.'^ , or the torus T'^ . We 
choose p > d' and a G ((p + d')/2p, 1). We denote X" = /}((— Aa?)") the fractional 
power space^ associated with the Laplacian operator A^v on hP{fi) with homogeneous 
Neumann boundary conditions. It is well-known^ that X" is continuously embedded in 
the Sobolev space W'^^'^^fl) and thus it is compactly embedded in C^{Q). Let r > 1 and 
/ G C(r2 X M X M'^',R). We consider the parabohc equation 

u{x, t) = Au{x, t) + f{x, u{x, t), Vu{x, t)) {x, t) eQ X (0, +oo) 

^{x,t) = {x,t) eOnx {0,+oo) (1.2) 

u(a;,0) = uo{x) G X° 

where u(t) is the time-derivative of u(t). 

Eq. (11. 2p defines a local dynamical system Sf{t) on X°' (see [IQ]) by setting Sf{t)uo = u{t). 

We assume moreover that there exist c G C°(]R+, M+), e > and k > such that / satisfies 

V/?>0, V^gR"^', s_up |/(x,^,OI<c(i?)(l + |e|'-^) 

and \/z G M, Vx G H, \z\> n => zf{x, z,0) <0 . 

Then, Eq. (II. 2p defines a global dynamical system in X'^ which admits a compact global 
attractor A (see [72]). 

The reader, which is not familiar with partial differential equations, may neglect all the 
technicalities about X'^, the Sobolev spaces and the parabolic equations in a first reading. 
The most important point is that Sf{t) is a dynamical system defined on an infinite- 
dimensional function space. Compared with the finite-dimensional case, new difficulties 
arise. For example, the existence of a compact global attractor requires compactness 
properties, coming here from the smoothing effect of (II. 2p . We also mention that, even if 
the backward uniqueness property holds, backward trajectories do not exist in general for 
(II. 2p . The reader interested in the dynamics of (II. 2p may consult [21|, [10], [SI], [Z2] or 

The purpose of this paper is to emphasize the different relationships between the dy- 
namics of (II. ip and (II. 2p . The correspondence is surprisingly perfect. It can be summarized 
by Table [H This correspondence has already been noticed for some of the properties of the 
table. We complete here the correspondence for all the known properties of the dynamics 
of the parabolic equation. Table [T] will be discussed in more details in Section [2] and, 
for cooperative systems, in Section HI Some of the properties presented in the table con- 
cerning finite-dimensional dynamical systems are trivial, other ones are now well-known. 



ODE 




PDE 


d= 1 

(Or more generally 

tridiagonal 

cooperative 

system of ODEs) 


• Gradient dynamics 

• Convergence to an equilibrium point 

• Automatic transversality of stable 
and unstable manifolds 

• Genericity of Morse-Smale property 

• Knowledge of the equilibrium points 
implies knowledge of the whole dynamics 

• Dimension of the attractor equal to the 
largest dimension of the unstable manifolds 

• Reahsation of the ODE in the PDE 


f^ = (0,l) 


d=2 
General case 

(Or more generally 
cyclic tridiagonal 

cooperative 
system of ODEs) 


• Poincare-Bendixson property 

• Automatic transversality of stable 
and unstable manifolds of two orbits 

if one of them is a hyperbolic periodic orbit 
or if both are equilibrium points 
with different Morse indices. 

• Non-existence of homoclinic orbits 
for periodic orbits 

• Genericity of Morse-Smale property 

• Reahsation of the ODE in the PDE 


= Ti 
General case 


d=2 
g radially 
symmetric 


• Automatic transversality of stable 
and unstable manifolds of 
equilibrium points and periodic orbits 

• No homochnic orbit 

• Knowledge of the equilibrium points 
and of the periodic orbits 

implies knowledge of the whole dynamics 

• Genericity of the Morse-Smale property. 

• Dimension of the attractor equal to the 
largest dimension of the unstable manifolds 

• Reahsation of the ODE in the PDE 


= Ti 
f{x,u,Vu) = f{u,Vu) 


d>3 


• Existence of persistent chaotic dynamics 

• Genericity of Kupka-Smale property 

• Reahsation of the ODE in the PDE 


dim(rj)> 2 


Any d 
g^VG 


• Gradient dynamics 

• Genericity of the Morse-Smale property 

• Realisation of a generic ODE in the PDE 


Any n 

f{x,u,Vu) = f{x,u) 



Table 1: the correspondence between the dynamics of vector fields and the ones of parabolic 
equations. 



However, the corresponding results for the parabohc equation are more involved and some 
of them are very recent. These properties are mainly based on Sturm-Liouville arguments 
and unique continuation properties for the parabolic equations as explained in Section [31 
The study of the dynamics generated by vector fields in dimension d > 3 is still a subject 
of research. Taking into account the correspondence presented in Table [1] should give a 
guideline for research on the dynamics of the parabolic equations. Some examples of open 
questions are given in Section O 

We underline that we only consider the dynamics on the compact global attractors. 
Hence, we deal with dynamical systems on compact sets. It is important to be aware of 
the fact that, even if the dimension of the compact global attractor A of the parabolic 
equation (11. 2p is finite, it can be made as large as wanted by choosing a suitable function 
/. This is true even if Q is one-dimensional. Therefore, all the possible properties of the 
dynamics of f ll.2p do not come from the low dimension of A but from properties, which 
are very particular to the flow of the parabolic equations. 

Finally, we remark that most of the results described here also hold in more general 
frames than (11. ip and (II. 2p . For example, W^ could be replaced by a compact orientable 
manifold without boundary. We could also choose for (II. 2p more general boundary condi- 
tions than Neumann ones, or less restrictive growing conditions for /. The domain Q may 
be replaced by a bounded smooth manifold. Finally, notice that the case Q = T"^ can be 
seen as ^2 = (0, l)'^ with periodic boundary conditions. 

2 Details and comments about the correspondence ta- 
ble 

We expect the reader to be familiar with the basic notions of the theory of dynamical 
systems and flows. Some definitions are briefly recalled in the glossary at the end of this 
paper. For more precisions, we refer for example to [IH], [SS], [SI], [20] or [M] for finite- 
dimensional dynamics and to [36], [IQ], [81] or [35] for the infinite-dimensional ones. 

We first would like to give short comments and motivations concerning the properties 
appearing in Table [H Notice that we do not deal in this section with the cooperative 
systems of ODEs. The properties of these systems are discussed in Section HI 

A generic property of the dynamics is a property satisfied by a countable intersection 
of open dense subsets of the considered class of dynamical systems. Generic dynamics 
represent the typical behaviour of a class of dynamical systems. For finite-dimensional 
flows, we mainly consider classes of the form (Tg(t))g^c^(^d^^dy The parameter is the vector 
field g, which belongs to the space C^(M'^,M'^) endowed with either the classical C^ or 
the C^ Whitney topology. Notice that the question wether or not a property is generic 
for g G C^(]R'^,R'^) for some r > 2 may be much more difficult than C^ genericity. We 
will not discuss this problem here. In some cases, we restrict the class of vector fields to 



subspaces of C^(M'^,M'^) such as radially symmetric, gradient vector fields or cooperative 
systems. In a similar way, for infinite-dimensional dynamics, we consider families of the 
type iSf(t))j:^(^i/Q^^y^^d'm, where C^{Q x R x M'^ ,M) is endowed with either the classical 
C^ or the C^ Whitney topology. For some results, we restrict the class of nonlinearities / 
to homogeneous ones or to ones, which are independent of the last variable ^. 

Poincare-Bendixson property and the convergence to an equilibrium or a periodic orbit 
are properties related to the following question: how simple are the a— and cj— limit sets 
of the trajectories ? For vector fields, the restriction of the complexity of the limit sets may 
come from the restriction of freedom due to the low dimension of the flow. As said above, 
there is no restriction on the dimension of the global attractor for the parabolic equations. 
The possible restrictions of the complexity of the limit sets come from particular properties 
of the parabohc equations, see Section |3J 

Hyperbolicity of equilibria and periodic orbits, tranversality of stable and unstable man- 
ifolds, Kupka-Smale and Morse-Smale properties are properties related to the question 
of stability of the local and global dynamics respectively. Morse-Smale property is the 
strongest one. It implies the structural stability of the global dynamics: if the dynamical 
system Tg{t) satisfies the Morse-Smale property, then for g close enough to g, the dynamics 
of Tg{t), restricted to its attractor Ag, are qualitatively the same as the ones of Tg(t) on 
Ag, see |60], [62] and |6T]. The same structural stability result holds for parabolic equa- 
tions satisfying the Morse-Smale property, see [36] , [35] and [59] . It is natural to wonder if 
almost all the dynamics satisfy these properties, that is if these properties are generic. 

The fact that the knowledge of the equilibria and the periodic orbits implies the knowl- 
edge of the whole dynamics may be studied at different levels. Two equilibria or periodic 
orbits being given, can we know if they are connected or not by a heteroclinic orbit ? Are 
two dynamics with the same equilibria and periodic orbits equivalent ? Is there a simple 
algorithm to determine the global dynamics from the position of the equilibria and the 
periodic orbits ? These questions are among the rare dynamical questions coming from 
the study of partial differential equations and not from the study of vector fields. Indeed, 
for finite-dimensional dynamical systems, the answers, either positive or negative, are too 
simple. In contrast, such kinds of results are probably among the most amazing ones for 
the dynamics of the parabolic equations. 

The persistent chaotic dynamics and the fact that the dimension of the attractor is equal 
to the largest dimension of the unstable manifolds, are related to the following question: 
how complicated may be the dynamics ? In general, the dimension of the attractor of a 
dynamical system may be larger than the largest dimension of the unstable manifolds. The 
classes of systems, where these dimensions automatically coincide, are strongly constrained, 
which in some sense implies a simple behaviour. On the contrary, chaotic dynamics have 
very complicated behaviour. Chaotic dynamics may occur through several phenomena, 
and the notion of chaotic behaviour depends on the authors. In this paper, "persistent 
chaotic dynamics" refers to the presence of a tranversal homoclinic orbit generating a Smale 



horseshoe (see [HB])- The persistent chaotic dynamics provide comphcated dynamics, which 
cannot be removed by small perturbations of the system. Such an open set of chaotic 
dynamics is a counter-example to the genericity of the Morse-Smale systems. 

The question of the realization of vector fields in the parabolic equations is as follows: 
a vector field g G C(]R'^, M"^) being given, can we find a function / and an invariant 
manifold M C L^(f2) such that the dynamics of the parabolic equation (11.21) restricted 
to M is equivalent to the dynamics generated by the vector field g 1 A positive answer 
to this question implies that the dynamics of the considered class of parabolic equations 
is at least as complicated as the dynamics of the considered class of vector fields. Such 
a realization result is very interesting since, on the opposite, the other properties stated 
in Table [1] roughly say that the dynamics of the parabolic equation (11.21) cannot be much 
more complicated than the ones of the corresponding class of finite-dimensional flows. One 
has to keep in mind that the manifold M, on which the finite-dimensional dynamics are 
realized, is not necessarily stable with respect to the dynamics of the parabolic equation. 
Typically, M cannot be stable if the finite-dimensional system contains a stable periodic 
orbit, since all periodic orbits of (II. 2p are unstable (see for example [S]). 

Now, we give short comments and references for the correspondences stated in Table [H 
• d = l and Vt = (0, 1) 

The dynamics generated by a one-dimensional vector field is very simple. Its attractor 
consists in equilibrium points and heteroclinic orbits connecting two of them. The exis- 
tence of these heteroclinic orbits is easily deduced from the positions of the equilibrium 
points. Moreover, these heteroclinic connections are trivially transversal. Finally, (II. ip 
is clearly a gradient system with associated Lyapounov functional F{y) = — J^ g{s)ds. 
As a consequence, the Morse-Smale property is equivalent to the hyperbolicity of all the 
equilibrium points, which holds for a generic one-dimensional vector field. 

The dynamics Sf{t) generated by (II. 2p for Q = (0, 1) is richer since its attractor may 
have a very large dimension. However, these dynamics satisfy similar properties. These 
similarities are mainly due to the constraints coming from the non- increase of the number of 
zeros of solutions of the linear parabolic equation (see Theorem 13. ip . Zelenyak has proved 
in [^ that Sf{t) admits an explicit Lyapounov function and thus that it is gradient. He 
also showed that the cu— limit sets of the trajectories consist in single equilibrium points. 
In Proposition 13.21 we give a short proof of this result, due to Matano. The fact that the 
stable and unstable manifolds of equilibrium points always intersect transversally comes 
from Theorem 13.11 and the standard Sturm-Liouville theory. This property has been first 
proved by Henry in [3S] and later by Angenent [2] in the weaker case of hyperbolic equilibria. 
As a consequence of the previous results, the Morse-Smale property is equivalent to the 
hyperbolicity of the equilibrium points and is satisfied by the parabolic equation on (0, 1) 
generically with respect to /. The most surprising result concerning (ll.2p on i7 = (0, 1) 
is the following one. Assuming that every equilibrium point is hyperbolic and that the 
equilibrium points ei,...,ep are known, one can say if two given equilibria e, and Cj are 



connected or not by a heteroclinic orbit. This property has been proved by Brunovsky 
and Fiedler in [12] for / = f{u) and by Fiedler and Rocha in [20] in the general case. The 
description of the heteroclinic connections is obtained from the Sturm permutation which 
is a permutation generated by the respective positions of the values ej(0) and ej(l) of the 
equilibrium points at the endpoints oi Q = (0, 1). The importance of Sturm permutation 
has been first underlined by Fusco and Rocha in [27] . We also refer to the work of Wolfrum 
[05] , which presents a very nice formalism for this property. Fiedler and Rocha showed 
in [22] that the Sturm permutation characterizes the global dynamics of (11.21) on (0,1). 
They proved in [2T] that it is possible to give the exact list of all the permutations which 
are Sturm permutations for some nonlinearity / and thus to give the list of all the possible 
dynamics of the parabolic equation on (0, 1). The fact that the dimension of the attractor 
is equal to the largest dimension of the unstable manifolds has been shown by Rocha in 
[52] . The previous works of Jolly [13] and Brunovsky [TT] deal with the particular case 
/ = f{u), but show a stronger result: the attractor can be embedded in a C^ invariant 
graph of dimension equal to the largest dimension of the unstable manifolds. Finally, let 
us mention that it is easy to realize any one- dimensional flow in an invariant manifold of 
the semi-flow generated by the one- dimensional parabolic equation. For example, in the 
simplest case of Neumann boundary conditions as in (11.21) . one can realize the flow of any 
vector field g as the restriction of the dynamics of the equation u = Au + g{u) to the 
subspace of spatially constant functions. 

• d = 2 and fi = T\ general case 

Even if they are richer than in the one- dimensional case, the flows generated by vector 
fields on M^ are constrained by the Poincare-Bendixson property (see the original works 
of Poincare [65j and Bendixson [7| or any textbook on ordinary differential equations). 
In particular, this constraint precludes the existence of non-trivial non-wandering points 
in Kupka-Smale dynamics. Due to the low dimension of the dynamics, the stable and 
unstable manifolds of hyperbolic equilibria or periodic orbits always intersect transversally 
if either one of the manifold corresponds to a periodic orbit or if the invariant manifolds 
correspond to two equilibrium points with different Morse indices. Moreover, there is no 
homoclinic trajectory for periodic orbits. Using these particular properties, Peixoto proved 
in [63j that the Morse-Smale property holds for a generic two-dimensional vector field. 

The first correspondence between two-dimensional flows and the dynamics of the para- 
bolic equation (II. 2p on the circle fi = T^ has been obtained by Fiedler and Mallet-Paret in 
[T9] . They proved that the Poincare-Bendixson property holds for (II. 2p on T^, by using the 
properties of the zero number (see Theorem 13. ip . The realization of any two-dimensional 
flow in a two-dimensional invariant manifold of the parabolic equation on the circle has 
been proved by Sandstede and Fiedler in [85|. Very recently, Czaja and Rocha have shown 
in [T7] that the stable and unstable manifolds of two hyperbolic periodic orbits always 
intersect transversally and that there is no homoclinic connection for a periodic orbit. The 
other automatic transversality results and the proof of the genericity of the Morse-Smale 



property have been completed by the authors in fl7| and 

• d = 2 and i7 = T^, radial symmetry and T^— equivariance 

When the vector field g satisfies a radial symmetry, the dynamics of the two-dimensional 
flow generated by (II. ip becomes roughly one-dimensional. The closed orbits consist in 0, 
circles of equilibrium points and periodic orbits being circles described with a constant 
rotating speed. The dynamics are so constrained that the closed orbits being given, it is 
possible to describe all the heteroclinic connections. Notice that no homoclinic connection 
is possible. We also underline that the Morse-Smale property is generic in the class of 
radially symmetric vector fields. 

If the two-dimensional radial vector fields are too simple to attract much attention, 
the corresponding case for the parabolic equation (11.21) on $1 = T^ with homogeneous 
nonlinearity /(x, m, dxu) = f{u, dxu) has been extensively studied. Since Theorem 13 . 1 I holds 
for (11. 2p with any one-dimensional domain fi, it is natural to expect results for (II. 2p on 
i7 = T^ similar to the ones obtained for (II. 2p on i7 = (0, 1). In particular, one may wonder if 
it is possible to describe the global dynamics of (II. 2p knowing the equilibria and the periodic 
orbits only. However, this property is still open for general non-linearities f{x,u,dxu) in 
the case ^2 = T^. Moreover, if one believes in the correspondence stated in this paper, 
one can claim that it is in fact false for a general nonlinearity f{x,u,dxu). Therefore, it 
was natural to first study the simpler case of homogeneous nonlinearities / = f{u,dxu). 
Indeed, the dynamics in this case are much simpler, in particular the closed orbits are either 
homogeneous equilibrium points e{x) = e G M, or circles of non-homogeneous equilibrium 
points, or periodic orbits consisting in rotating waves u{x,t) = u{x — ct) (notice the 
correspondence with the closed orbits of a radially symmetric two-dimensional fiow). This 
property is a consequence of the zero number property of Theorem 13.11 and has been proved 
in [1] by Angenent and Fiedler. The works of Matano and Nakamura [SI] and of Fiedler, 
Rocha and Wolfrum [23] show that the unstable and stable manifolds of the equilibria and 
the periodic orbits always intersect transversally and that no homoclinic orbit can occur. 
Moreover, in [23] , the authors give an algorithm for determining the global dynamics of the 
parabolic equation (II. 2p on fi = T^ with homogeneous nonlinearity / = f{u,dxu). This 
algorithm uses the knowledge of the equilibria and the periodic orbits only. In [83] , Rocha 
also characterized all the dynamics, which may occur. Due to the automatic transverse 
intersection of the stable and unstable manifolds and due to the possibility of transforming 
any circle of equilibrium points into a rotating periodic orbit (see [23]), one can show that 
the Morse-Smale property holds for the parabolic equation on T^ for a generic homogeneous 
nonlinearity f{u,dxu) (see j47j). Finally, the realization of any radially symmetric two- 
dimensional fiow in the dynamics of (II. 2p on T^ for some / = f{u, d^u) and the fact that 
the dimension of the attractor is equal to the largest dimension of the unstable manifolds 
are shown in [35] . 

• d>?, and dim(fi)> 2 

The genericity of the Kupka-Smale property for vector fields in M*^, c? > 3, has been proved 



independently by Kupka in [50j and by Smale in |S7] . Their proofs have been simphfied by 
Peixoto in [Bl] (see [T] and [61J). The strong difference with the lower dimensional vector 
fields is that, when d > 3, (11. ip may admit transversal homoclinic orbits consisting in the 
transversal intersection of the stable and unstable manifolds of a hyperbolic periodic orbit. 
The existence of such an intersection is stable under small perturbations and yields a Smale 
horseshoe containing an infinite number of periodic orbits and chaotic dynamics equivalent 
to the dynamics of the shift operator, see [SS]. Therefore, the Morse-Smale property cannot 
be dense in vector fields on M*^ with d > 3. Even worse, the set of vector fields, whose 
dynamics are structurally stable under small perturbations, is not dense (notice that this 
set contains the vector fields satisfying the Morse-Smale property). Indeed, as shown in 
[32] . there exist sets of vector fields gx, depending smoothly of a real parameter A, such that, 
for each A, Tg^{t) admits a Lorenz attractor, whose dynamics are qualitatively different for 
each value of A. The possible presence of other chaotic dynamics such as Anosov systems 
or wild dynamics is also noteworthy, see |5], [55] and [9]. For the interested reader, we 
refer to [^ or [89]. 

The genericity of the Kupka-Smale property for the parabolic equation (ll.2p on any 
domain Q is proved by Brunovsky and both authors in \IS\. There exist several results 
concerning the embedding of the finite-dimensional flows into the parabolic equations. 
Polacik has shown in [68] that any ordinary differential equation on M*^ can be embedded 
into the flow of (ll.2p for some / and for some domain fi C M*^. The constraint that 
the dimension of Q is equal to the dimension of the imbedded flow is removed in [69] . 
however the result concerns a dense set of flows only. A similar result has been obtained 
by Dancer and Polacik in [18j for homogeneous nonlinearities f{u, Vm) (see also [73]). These 
realization results imply the possible existence of persistent chaotic dynamics in the flow 
of the parabolic equation (II. 2p as soon as Q has a dimension larger than one: transversal 
homoclinic orbits, Anosov flows on invariant manifolds of any dimension, Lorenz attractors 
etc... 

• Gradient case 

When g is a gradient vector field VG with G G C^(]R'^,R), then —G is a strict Lyapounov 
function and (11.10 is a gradient system. In this case, the Kupka-Smale property is equiva- 
lent to the Morse-Smale property. The genericity of the Morse-Smale property for gradient 
vector fields has been obtained by Smale in [86] . 

In the case where the nonlinearity / e C^{Q x R, R) (that is, / = f{x,u) does not 
depend on Vm), the parabohc equation (II. 2p admits a strict Lyapounov function given 
by E{u) = J^(^^\Vu{x)\'^ — F{x,u{x))) dx, where F{x,u) = J^ f{x,s)ds is a primitive 
of /, and hence generates a gradient system. Brunovsky and Polacik have shown in [H] 
that the Morse-Smale property holds for the parabolic equation, generically with respect 
to f{x,u). It is noteworthy that the Morse-Smale property is no longer generic if one 
restricts the nonlinearities to the class of homogeneous functions / = f{u) (see [71]). 
Polacik has shown in [7D] that any generic gradient vector field of M'^ can be realized in 
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the flow of tlie parabolic equation f ll.2p on a bounded domain of M^ witli an appropriate 
nonlinearity f{x,u). Tlie paper [TD] also contains the realization of particular dynamics 
such as non-transversal intersections of stable and unstable manifolds. 

• Caveat: general ODEs or cooperative systems? 

In Table 1, we have given the striking correspondence between the flow generated by Eq. 
(11. ip and the semifiow generated by the parabolic equation (11.21) . In addition, we have 
pointed out that some classes of cooperative systems are also involved in this correspon- 
dence. In fact, the reader should be aware that the dynamics of (II. 2p is much closer to the 
ones of a cooperative system than to the ones of the general vector field (II. ip . Indeed, the 
semifiow Sf{t) generated by the parabolic equation (II. 2p belongs to the class of strongly 
monotone semiflows, which means that this semifiow has more constraints than the flow 
Tg{t) generated by a general vector field g (see Section H]). That is why, it could be more 
relevant to write Table 1 in terms of cooperative systems only (for example, by replacing 
the case of the general ODE with d > 3 by the case of a cooperative system of ODEs in 
dimension d > 4). However, we have chosen to mainly write Table 1 in terms of general 
ODEs for several reasons: 

- as far as the properties stated in Table 1 are concerned, there is no difference between 
the dynamics of a general ODE and the ones of a parabolic PDE, 

- the dynamics of general ODEs are common knowledge, whereas speaking in terms of 
cooperative systems may not give a good insight of the dynamics of (II. 2p . 

- not all the properties stated in Table 1 are known for the class of cooperative systems 
(for example the genericity of Kupka-Smale property is not yet known for li > 4). 

3 Zero number and unique continuation properties 
for the scalar parabolic equation 

The results presented in Table [Hand in Section [2] strongly rely on properties specific to the 
parabolic equations. The purpose of this section is to give a first insight of these particular 
properties and of their use, to the reader. 

Dynamical systems generated by vectors fields are flows on Mf^, whereas the phase-space 
of the parabolic equation is an infinite-dimensional space X". It is important to be aware 
of the fact that the parabolic equations generate only a small part of all possible dynamical 
systems on the Banach space X". On one hand, this implies less freedom in perturbing 
the dynamics and hence in obtaining density results. In particular, whereas one can easily 
construct perturbations of a vector field g which are localized in the phase space W^, the 
perturbations of the nonlinearity / act in a non local way on X"' (many different functions 
u can have the same values of u and Vtt at a given point x). Therefore, it is important 
to obtain unique continuation results in order to find values (x, m, Vm), which are reached 
only once by a given periodic, heteroclinic or homoclinic orbit. On the other hand, the 
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small class of dynamics generated by the parabolic equations admits special properties. 
These properties may in particular yield the constraints, which make the dynamics similar 
to the ones of low- dimensional vector fields. 

The scalar parabolic equation in space dimension one [Q = (0, 1) or T^) satisfies a 
very strong property: the number of zeros of the solutions of the linearized equation is 
nonincreasing in time. This property is often called Sturm property since its idea goes 
back to Sturm [94j in 1836. There are different versions of this result, which have been 
proved by Nickel [SZl, Matano [521 E3], Angenent and Fiedler [31 S] (see also [2S] for a 
survey). By similar technics, a geometrical result on braids formed by solutions of the 
one-dimensional parabolic equation is obtained in [3T] . 



Theorem 3.1. Let Q = (0, 1) with Neumann boundary conditions or Q = T^ . Let T > 0, 

a G W^'^iU X [0,T],M) and b G L~(n x [0,T],M). Letv -.Ux (0,T) ^ R be a bounded 
non-trivial classical solution of 

dtv = dl^v + a{x,t)d:,v + b{x,t)v , (x, t) G fi x (0,T) . 

Then, for any t G (0, T), the number of zeros of the function x G fi i— )■ f (x, t) is finite and 
non-increasing in time. Moreover, it strictly decreases at t = tQ if and only if x \-^ f (x, to) 
has a multiple zero. 

Theorem 13.11 is the fundamental ingredient of almost all the results given in Table [1] in 
the cases Vt = (0, 1) and Q = T^. It can be used either as a strong comparison principle or as 
a strong unique continuation property, as shown in the following examples of applications. 
General surveys can be found in [21], [M] and [55] . 

In the first application presented here. Theorem 13.11 is used as a strong maximum 
principle. In some sense, it yields an order on the phase space which is preserved by the 
flow. This illustrates how Theorem 13.11 may imply constraints similar to the ones of low- 
dimensional vector fields. The following result was first proved in [HS] and the proof given 
here comes from [52] (see also |24j). 

Proposition 3.2. Let Q = (0,1), let uq G X" and let u{x,t) be the corresponding solu- 
tion of the parabolic equation (11.21) with homogeneous Neumann boundary conditions. The 
u— limit set of uq consists of a single equilibrium point. 

Proof: We first notice that v{x,t) = dtu{x,t) satisfies the equation 

dtv{x, t) = dl^v{x, t) + f^{x, u{x, t), dxu{x, t))v{x, t) + /a^„(x, u{x, t), d^u{x, t))d^v{x, t) . 

Due to the Neumann boundary conditions, we have ^^^^(O,^) = dl^u{0,t) = dxv{0,t) = 
for all t > 0. In particular, as soon as v{0, t) = 0, v{t) has a double zero at x = 0. Due to 
Theorem 13. 1[ either f is a trivial solution, that is f = for all t, and u is an equilibrium 
point, or v{0,t) vanishes at most a finite number of times since v{t) can have a multiple 
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zero only a finite number of times. Assume tliat u is not an equilibrium, then u{0,t) must 
be monotone for large times and thus converges to a G M. Any trajectory w in the w— limit 
set of uo must hence satisfy w{0,t) = a for all t. Therefore, dtw{0,t) = for all t and 
dtw{0,t) has a multiple zero at x = for all times. Using Theorem 13. ![ we deduce as 
above that w is an equilibrium point of (11. 2p . But there exists at most one equilibrium 
w satisfying w{0) = a and the Neumann boundary condition dxw{0) = 0. Therefore, the 
w— limit set of Uq is a single equilibrium point w. D 



The second application comes from |17] . It shows how Theorem 13.11 can be used as a 
unique continuation property. This kind of property roughly says that, if two solutions co- 
incide too much near a point {xq, to), then they must be equal everywhere. The motivation 
beyond this example of application is the following. We consider a time-periodic solution 
of (11. 2p on r2 = T^. The problem is to find a perturbation of the nonlinearity /, which 
makes this periodic orbit hyperbolic. As enhanced above, such a perturbation is nonlocal 
in the phase space of (II. 2p . To be able to perform perturbation arguments, it is important 
to show that one can find a perturbation of / which acts only locally on the periodic orbit. 
To this end, one proves the following result. 

Proposition 3.3. Let p{x,t) be a periodic orbit of (II. 2p on Q = T^. Let T > be its 
minimal period. Then, the map 

{x,t) e T^ X [0,T) h-> {x,p{x,t),d,p{x,t)) 

is one to one. 

Proof: Assume that this map is not injective. Then there exist Xq, to ^ [0,^) ^-^id 
ti G [0,T), to 7^ ti such that 

p(a;o,to) =p(xo,ti) and d^p{xo,to) = d^p{xo,ti) . 

The function v{x, t) = p{x, t + ti — to) — p{x, t) is a solution of the equation 

dtv{x, t) = dl^v{x, t) + a(x, t)v{x, t) + 6(x, t)d^v{x, t) , 

where a{x,t) = f^ f'^{x,p{x,t) + s{p{x,t + ti — to) — p{x,t)),dxp{x,t + ti — to))c?s and 
b{x,t) = Jq fl^{x,p{x,t),dx{p{x,t) + s{p{x,t + ti — to)—p{x,t)))ds. Moreover, the function 
v{x,t) satisfies f(a:o,to) = and dxv{xo,to) = and does not vanish everywhere since 
1*^1 ~ ^o| < T. Due to Theorem 13.1 1 the number of zeros of v{t) drops strictly at t = to and 
never increases. However, v{t) is a periodic function of period T, and thus, its number of 
zeros is periodic. This leads to a contradiction and proves the proposition. D 
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In a domain Q of dimension d' > 2, there is no known counterpart for Theorem 13.11 
as shown in [2H]- In particular, Proposition 13.31 does no longer hold. However, to be 
able to construct relevant perturbations of periodic orbits, one needs a result similar to 
Proposition 13. 3[ even if weaker. The following result can be found in [T3] . Its proof is 
based on a generalization of the arguments of [37] and on unique continuations properties 
of the parabolic equations. 

Theorem 3.4. Let p{x,t) be a periodic orbit of (11. 2p with minimal period T > 0. There 
exists a generic set of points {xo,to) e fi x [0, T) such that if t E [0,T) satisfies p{xo,t) = 
p{xo,to) and Vp(xo,t) = Vp(xo,to), then t = to- 

4 Cooperative systems of ODEs 

We consider a system of differential equations 

yit) = givit)) , i/(0) = yoGK^, (4.1) 

where g = {gi)i=i...N is a C^ vector field. 

Due to the analogy with biological models, the following definitions are natural. We say 

that (14. ip is a cooperative (resp. competitive) system if for any y G M^ and % ^ j-, -^{y) 

is non-negative (resp. non-positive) and the matrix {^){y) is irreducible i.e. it is not a 
block diagonal matrix (the simpler assumption that all the coefficients -^{y) are positive 
is sometimes made instead of the irreducibility) . We say that (14. ip is a tridiagonal system 
if ^ = for I? — j| > 2 and a cyclic tridiagonal sysiera if the indices i and j are considered 

modulo N , i.e. if, in addition, we allow tj^ and -#^ to be non-zero. For the reader 

' ' ' aVN ayi 

interested in cooperative systems, we refer to |93] . 

In this section, we only consider cooperative systems. However, notice that, by changing 
t into —t or yi into — ^j, we obtain similar results for competitive systems and for systems 
with different sign conditions. 

The dynamics of cooperative systems may be as complicated as the dynamics of general 
vector fields. Indeed, Smale has shown in [90] that any vector field in M^^^ can be realized 
in a invariant manifold of a cooperative system in R^. Notice that this realization result 
implies that any one- dimensional vector field can be imbedded in a tridiagonal coopera- 
tive system and any two-dimensional vector field can be imbedded in a cyclic tridiagonal 
cooperative system. This explains why we present the tridiagonal cooperative systems in 
Table [1] as generalization of one- and two-dimensional vector fields. 

However, the dynamics of a cooperative system (14. ip is really different from the ones of 
the general ODE (II. ip since a cooperative system generates a strongly monotone fiow, that 
is, a fiow which preserves a partial order. It is noteworthy that the semifiow Sf{t) generated 
by the parabolic equation (11.21) also belongs to the class of strongly monotone semifiows (it 
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preserves the order of X" induced by the classical order of C^{Q)). Therefore, the semiflow 
of fll.2p is much closer to the one of the cooperative system (14. 11) . both admitting more 
constraints than the flow Tg{t) generated by a general vector field g. In [H] and [H] for 
example, Hirsch has shown that almost all bounded trajectories of a strongly monotone 
semiflow are quasiconvergent, that is, their w-limit sets consist only of equilibria. More 
precisely, all initial data, which have bounded nonquasiconvergent trajectories, form a 
meager subset (that is, the complement of a generic subset) of the phase space. Later, 
in |67], Polacik has proved that the set of all initial data uq G X", which have bounded 
nonconvergent trajectories in the semiflow of the parabolic equation (11.21) . is meager in 
X°. 

Moreover, since the works of Hirsch and Smillie, it is known that the dynamics of 
cooperative systems, which are in addition tridiagonal, are very constrained in any dimen- 
sion N. Indeed, in [12], |13] and [HI], strong properties of the limit sets of cooperative 
systems are proved. In particular, any three-dimensional cooperative system satisfies the 
Poincare-Bendixson property and the trajectory of any tridiagonal cooperative system con- 
verges to a single equilibrium point. Inspired by the articles of Henry and Angenent about 
the parabolic equation on (0,1), Fusco and Oliva (see [2S]) showed a theorem similar to 
Theorem 13.11 (see [H2] for a more general statement). 

Theorem 4.1. Let M he the set of vector y G M^ such that, for all i = 1 . . .N, either 
Hi j^ or Hi = and yi^iVi+i < (where y^ = Vn+i = ^)- For every y G N , we set N{y) to 
be the number of sign changes for yi, when i goes from 1 to N . Let y{t) ^ be a solution 

of 

y{t) = A{t)y{t) , (4.2) 

where A G C°(]R, A^Ar(]R)) satisfies Aij{t) > for allt eR and all i^ j. 

Then, the times t where y{t) ^ J\f are isolated and, if yit^) ^ J\f, then, for every e > 

small enough, N{y{t + e)) < N{y{t — e)). 

In other words, the number of sign changes of the solutions of the linear equation (14. 2 p 
is non- increasing in time and strictly drops at to if ^^^ only if y{to) has in some sense a 
multiple zero. The parallel with Theorem 13.11 is of course striking. Using Theorem 14. H 
Fusco and Oliva have shown that the stable and unstable manifolds of equilibrium points 
of a tridiagonal cooperative system always intersect transversally. As a consequence, the 
Morse-Smale property is generic in the class of tridiagonal cooperative systems. 

Theorem 14.11 also holds for cyclic tridiagonal cooperative systems, see [2S] and |92] . 
Using this fundamental property, Fusco and Oliva have shown in [20] that the stable and 
unstable manifolds of periodic orbits of cyclic tridiagonal cooperative systems always in- 
tersect transversally. In addition, Mallet-Paret and Smith have shown in [51] that cyclic 
tridiagonal cooperative systems satisfy the Poincare-Bendixson property. Notice that, fol- 
lowing [17] and [IB], one should be able to prove the genericity of the Morse-Smale property 
for cyclic tridiagonal cooperative systems. This has been proved very recently by Percie 
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du Sert (see 



Considering all these results, it is not surprising that there exists a parallel between 
tridiagonal cooperative systems and the parabolic equation on (0,1). Indeed, consider a 
solution V of the linear one-dimensional parabolic equation 

v{x,t) = dl^v{x,t) + a{x,t)d^v{x,t) + b{x,t)v{x,t) (a;, t) G M x IR+ (4.3) 

We discretize the segment (0, 1) by a sequence of points Xk = {k — 1)/{N — 1) with 
k = 1 . . . N. The natural approximation of v is given by yk ~ v{xk) solution of 

. ,,x yk+i{t) - 2yk{t) + yk-i{t) yk+i - yk , , (.. ,.. , . .. 

Vkit) = — + ak{t) + bk{t)yk{t) (4.4) 

where afc(t) = a{xk,t), bk(t) = b{xk,t) and h = 1/{N — 1). If h is small enough, (14.41) is a 
tridiagonal cooperative system. The relation between Theorems 13.11 and 14.11 is obvious in 
this framework. 

5 Beyond Kupka-Smale and other open problems 

One of the main goals of the study of dynamical systems is to understand the behaviour of a 
generic dynamical system. The most recent results concerning the parabolic equations are 
related to the genericity of Kupka-Smale property. However, such a property cannot give 
a good insight of the complex and chaotic dynamics that may be generated by homoclinic 
connections. For this reason, the study of finite-dimensional flows has been pursued further 
the Kupka-Smale property and is still in progress. The corresponding results should serve 
as a guideline for the study of the flow generated by the parabolic equation (II. 2p . 

As for vector fields, going beyond the Kupka-Smale genericity property for scalar 
parabolic equations will be an important and difficult problem. For vector fields, one 
of the main steps for this purpose is Pugh's closing lemma: if p is a non-wandering point 
of the dynamical system Tg[t) generated by (II. ip . then there exists a C^— perturbation g of 
g such that p is a periodic point of Tg{t) (the case of a C^— perturbation with r > 2 is still 
open). The proof of Pugh in [71] concerns discrete dynamical systems. It has been adapted 
to the case of fiows by Pugh and Robinson in [TB] (see also [70] for an introduction to the 
proof). A direct consequence of Pugh closing lemma is the general density theorem: for a 
generic finite-dimensional fiow, the non-wandering points are is the closure of the periodic 
points (see [75] and [80]). Other connecting lemmas have been proved by Hayashi [38] and 
Bonatti and Crovisier [S]. They enable a better understanding of generic dynamics. For 
example, the class of finite-dimensional dynamical systems which either satisfy the Morse- 
Smale property or admit a transversal homoclinic connection is generic (see [77], [10] and 
[16] for discrete dynamical systems in dimensions d = 2, d = 3 and d> 4 respectively, and 
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see [H] for three-dimensional flows). Obtaining similar results for the flow of the parabolic 
equation should be a very interesting and difficult challenge. 

Other interesting open problems concern the realization of flnite-dimensional dynamics 
in the semiflow of parabolic equations. Indeed, we only know that one can realize the 
dynamics of a dense set of general ODEs in the flow of a parabolic equation (11.21) on a two- 
dimensional domain. One may wonder if it is possible to realize the dynamics of all ODEs. 
Since the parallel between parabolic equations and cooperative systems is stronger, the 
following strong realization conjecture may be more plausible: any flow of a cooperative 
system of ODEs can be realized in an invariant manifold of the flow of a parabolic equation 
(11. 2p on a two-dimensional domain. 

Finally, the genericity of the Morse- and Kupka-Smale properties is also an interesting 
problem for other classes of partial differential equations. The genericity of the Morse- 
Smale property is known for the wave equations u + 'ju = Au + f{x,u) with constant 
damping 7 > (see [15]) and with variable damping 7(x) > in space dimension one (see 
[16]). We recall that, in both cases, the associated dynamical system is gradient. Nothing 
is known for other classes of PDEs, in particular for the equations of fluids dynamics and 
for systems of parabolic equations U = AU + f{x,U), with U{x,t) G M^. In all these 
cases, the main problem consists in understanding how the perturbations act on the phase 
plane of the PDE. Either one proves unique continuation results similar to Theorem 13.41 
in order to be able to use local perturbations of the flow (as in [17], [H] and [13]), or one 
uses particular non-local perturbations in a very careful way (as in [H], [15] and [16]). 

Glossary 

In this section, S{t) denotes a general continuous dynamical system on a Banach space X. 
An orbit of S{t) is denoted by x{t) = S(t)xo with t E I, where / = [0, +00), / = (—00, 0] 
or / = (—00, +00) in the case of a positive, negative or global trajectory respectively. 
Compact global attractor: if it exists, the compact global attractor A of S{t) is a 
compact invariant set which attracts all the bounded sets of X. Notice that A is then the 
set of all the bounded global trajectories. See [33] . 

a— and w— limit sets: let Xq G X. The a— limit set a{xQ) and the w— limit set a;(xo) of 
Xq are the sets of accumulation points of the negative and positive orbits coming from xq 
respectively. More precisely, 

a{xo) = {x E X / 3(t„)„gN, t„ y —00 and a negative trajectory x{t) 

such that x(0) = Xq and x(t„) ;■ x} 

u{xo) = {x G X / 3(trt)nGN, tn > +00 such that S(tn)xo > x} 

The limit sets a(xo) and w(xo) are non-empty connected compact sets. 
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Homoclinic or heteroclinic orbit: let x{t) = S{t)xQ be a global trajectory of S{t). 
Assume that the a— and w— limit sets of Xq exactly consists in one orbit, denoted a;_(t) 
and x^{t) respectively, this orbit being either an equilibrium point or a periodic orbit. The 
trajectory x{t) is said to be a homoclinic orbit if X-{t) = x+{t) and a heteroclinic orbit if 
x_(t)^x+(t). _ 

Backward uniqueness property: S{t) satisfies the backward uniqueness property if for 
any time tg > and any trajectories Xi{t) and X2{t), a;i(to) = 2^2(^0) implies Xi(t) = X2{t) 
for all t G [0,to]- Notice that this does not mean that S{t) admits negative trajectories. 
Hyperbolic equilibrium points or periodic orbits: an equilibrium point e of S{t) 
is hyperbolic if the linearized operator x h-)- DeS{l)x has no spectrum on the unit circle. 
Let p{t) be a periodic solution of S{t) with minimal period T. For each a; G X, we 
denote t 1— )■ n(t, 0)a; the corresponding trajectory of the linearization of S{t) along the 
periodic solution p{t). Then, p{t) is said hyperbolic if the linear map x ^-> n(T, 0)x has no 
spectrum on the unit circle except the eigenvalue 1 which is simple. Remark that then, 
for any integer k ^ 0, the linear map x 1— )• Il{kT, 0)x has no spectrum on the unit circle 
except the eigenvalue 1 (which is simple). 

Stable and unstable manifolds: let e be a hyperbolic equilibrium point of S{t). There 
exists a neighbourhood A/" of e such that the set 

^loci^) — {^0 G X , 3 a negative trajectory x{t) with x(0) = Xq 

and, Vt < 0, x{t) G M} 

is a submanifold of X, in which all negative trajectories converge to e when t goes to —00. 
The manifold Wl^^{e) is called the local unstable manifold of e. Pushing Wl^^{e) by the flow 
S(t), one can define the (global) unstable set W^{e) = Ut>oS(t)Wl^^{e), which consists in 
all the negative trajectories converging to e when t goes to — cxd. This unstable set W^{e) 
is an immersed submanifold under suitable properties. For instance, backward uniqueness 
properties are needed to extend the manifold structure. 
In the same way, one defines the local stable manifold 

Wfje) = {xo G X , Vt > 0, S{t)xo G A^} 

= {xo G X , Vt > 0, S{t)xo G TV and S{t)xo > e} . 

i— i>+oo 

General partial differential equations (and parabolic equations in particular) do not admit 
negative trajectories. Therefore, it is less easy to extend the local stable manifold to a 
global stable manifold. However, one can define the stable set W^{e) of e as follows 

W'{e) = {xo G X , S{t)xo > e} . 

i— >-+oo 

Under suitable additional properties (which are satisfied by the parabolic equation (II. 2p ). 
one can show that P^**(e) is also an immersed submanifold. For instance, backward unique- 
ness properties of the adjoint dynamical system S*{t) on X* and finite-codimensionality 
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of VF;^^(e) are needed (see [10] for more details). 

li p{t) is a hyperbolic periodic orbit, one defines its unstable and local stable manifolds in 

the same way. See for example [61] for more details. 

Non- wandering set: a point Xq G X is no n- wandering if for any neighbourhood M 3 Xq 

and any time to > 0, there exists t >tQ such that S{t)Af fl A/" 7^ 0. 

The Kupka-Smale and Morse-Smale properties: S{t) satisfies the Kukpa-Smale 

property if all its equilibrium points and periodic orbits are hyperbolic and if their stable 

and unstable manifolds intersect transversally. It satisfies the Morse-Smale property if in 

addition its non-wandering set consists only in a finite number of equilibrium points and 

periodic orbits. We refer to [6IJ for more precisions on these notions. 

Gradient dynamical systems: S{t) is gradient if it admits a Lyapounov functional, that 

is a function $ G C°(X, M) such that, for all Xq ^ X, t \-^ $(S'(t)xo) is non-increasing and 

is constant if and only if xq is an equilibrium point. We recall that a gradient dynamical 

system does not admit periodic or homoclinic orbits. 

Cooperative system of ODEs: see Section SJ 

Generic set and Baire space: a generic subset of a topological space F is a set which 

contains a countable intersection of dense open subsets of Y. A property is generic in Y if 

it is satisfied for a generic set of Y. The space Y is called a Baire space if any generic set 

is dense in Y. In particular a complete metric space is a Baire space. 

Whitney topology: let A; > and let M be a Banach manifold. The Whitney topology 

on C^{M,M.) is the topology generated by the neighbourhoods 

{geC\M,R) , \D'f{x)-D'g{x)\ <6{x), Vz G {0, 1, . . . , A;}, Vx e M} , 

where / is any function in C'^(M, M) and 6 is any positive function in C^(M, (0, +00)). 
Notice that C^{M, M) endowed with the Whitney topology is a Baire space even if it is not 
a metric space when M is not compact. We refer for instance to |30] . 
The fractional power space X": let A be a positive self-adjoint operator with compact 
inverse on L^(fi). Let (A^,) be the sequence of its eigenvalues, which are positive, and let 
(ifn) be the corresponding sequence of eigenfunctions, which is an orthonormal basis of 
lJ{Q). For each a G M, we define the fractional power of A by A(^^ Cnfn) = J2n (^n^nfn- 
In particular, A^ = Id and A^ = A. The space X" is the domain of A" that is 

X" = {ipe V{n) , if = ^c^y^n such that (c„A^) G l^N)} . 

n 

It is possible to define fractional powers of more general operators, called sectorial opera- 
tors, see [lO] . 

The Sobolev space W'^''^{fl): if s is a positive integer, iy*'^(f2) is the space of (classes 
of) functions / G hP{Q), which are s times differentiable in the sense of distributions and 
whose derivatives up to order s belong to L^(r2). It is possible to extend this notion to 
positive numbers s which are not integers by using interpolation theory. 
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Unique continuation properties: let us consider a partial differential equation on Q 
and let u{x,t) be any solution of it. A unique continuation property for this PDE is a 
result stating that if u{x,t) vanishes on a subset of 17 x R_|_ which is too large in some 
sense, then u{x,t) must vanish for all {x,t) in i7 x ]R+. 
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